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A novel route to a finite center-of-mass momentum pairing state; current driven 

FFLO state 

Hyeonjin Dohy Matthew Song, and Hae- Young KeqJ 
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(Dated: February 6, 2008) 

The previously studied Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state is stabilized by a magnetic 
field via the Zeeman coupling in spin-singlet superconductors. Here we suggest a novel route to 
achieve non-zero center-of-mass momentum pairing states in superconductors with Fermi surface 
nesting. We investigate two-dimensional superconductors under a uniform external current, which 
leads to a finite pair-momentum of q^. We find that an FFLO state with a spontaneous pair- 
momentum of q^ is stabilized above a certain critical current which depends on the direction of the 
external current. A finite qs arises in order to make the total pair-momentum of qt(= qs -I- qe) 
perpendicular to the nesting vector, which lowers the free energy of the FFLO state, as compared 
to the superconducting and normal states. We also suggest experimental signatures of the FFLO 
state. 

PACS numbers: 74.25.Dw,74.25.Sv,74.81.-g 



Introduction: Fulde and Ferrell|l|, and Larkin 
and Ovchinnikov 2] predicted that Cooper pairs with 
a nonzero center-of-mass momentum can be stabilized 
under a finite magnetic field via the Zeeman couphng, 
when the Zeeman term dominates over the orbital ef- 
fect. Recently, there has been growing interest in the- 
oretical and experimental studies of the Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) state in spin-singlet super- 
conductors and trapped cold fermionic atoms with un- 
equal densities [3J. It has been suggested theoretically 
that the organic superconductors like k-(ET)2 salts and 
A-(ET)2 salts, heavy-fermion superconductors, and high 
Tc cuprates are promising candidates for observing the 
FFLO state. The formation of a possible FFLO state 
has been very recently inferred from specific heatl4| and 
magnetization measurements 13 on one of heavy fermion 
systems, CeCoIus known to be a d-wave superconductor. 

The main idea of the previously studied FFLO state is 
based on the well known BCS theory that Cooper pairs in 
spin-singlet superconductors are made of fcrmions with 
opposite spins. When the magnetic field is acting on the 
spin of the electrons via the Zeeman coupling, electrons 
tend to polarize along the direction of the magnetic field 
to gain polarization energy, while the pairing of opposite 
spins is favorable for condensation energy. As a result of 
the competition, the superconducting state undergoes a 
transition to the FFLO state, and the system eventually 
enters the normal state, as the magnetic field is further 
increased. It has been shown that the FFLO state can 
lead to an enhancement of the critical magnetic field of 
up to 2.5 times the Pauli paramagnetic limit. It has also 
been suggested that the FFLO state can be realized when 
two species of fermion have different densities in trapped 
cold fermionic atoms. I3I 
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In this paper, we suggest a novel route to achieve the 
FFLO state. We study two dimensional superconductors 
in the presence of uniform external current, jg. When a 



FIG. 1: (Color online) The thick arrow indicates a nesting 
vector of a Fermi surface on a square lattice. When an exter- 
nal current with a momentum qe is applied along the a;-axis, 
a spontaneous pair-momentum of qa is generated in such a 
way that the total pair-momentum, qt(= qe + qs) is perpen- 
dicular to the nesting vector. The region of Cooper pairing 
between electrons with momentum k -I- qt/2 and — k -I- qt/2 
is indicated by the dark bars . This pairing on the Fermi sur- 
face has the dominant contribution to lowering the free energy 
of the FFLO state, as compared to the superconducting and 
normal states. 



current fiows. Cooper pairs acquire a finite external pair- 
momentum of Qe, and the superconducting state becomes 
unstable towards the normal state at a critical current. If 
the underlying Fermi surface has nesting, we find that the 
strength of this critical current depends on the direction 
of the current. In the case of s-wave superconductors, the 
critical current is maximal when the current is perpendic- 
ular to the nesting vector. We show that the FFLO state 
with a spontaneous pair-momentum of q^ is stabilized, 
when the current is not perpendicular to the nesting vec- 
tors. For example, a spontaneous pair-momentum, q^ is 



induced perpendicular to the current, when the current is 
along the x-axis of the square lattice as shown in Fig. ^ 
The total pair- momentum of qj (= qe + Qs) is determined 
to maximize the Cooper pairing between electrons with 
momentum k + qt/2 and — k + qt/2 (in a lab frame). 
This pairing on the Fermi surface has the dominant con- 
tribution to lowering the free energy of the FFLO state, 
as compared to the superconducting and normal states. 
Superconductors in the presence of a uniform 
supercurrent: We consider a sample with thickness of 
d <C ^, where ^ is the coherence length. Under this condi- 
tion, the magnitude of the superconducting gap and the 
current are uniform across the sample. When a uniform 
current flows, the Cooper-pair acquires a center-of-mass 
momentum of qe. For an infinitesimal current there will 
be no contribution from the quasiparticles, and the su- 
percurrent will be simply proportional to q^- However, 
the change of the quasiparticle spectrum in the presence 
of a supercurrent known as a Doppler shift, is written as 
follows. 




FIG. 2: (Color online) The direction dependence of the crit- 
ical current value (in arbitrary unit) is shown as the thick 
dashed line. The thin solid line is the Fermi surface of the 
tight binding model on the square lattice. The critical current 
has its maximum, j'c'^^ when the current is perpendicular to 
the nesting vector, while its minimum is j™'" along the a; (or 
2/)-direction. 



^k = ^1^ + Vk • q J2 



(1) 



where E^ is the quasiparticle spectrum in the absence 
of the current, and Vk is the quasiparticle velocity. The 
total supercurrent as a function of qe rises linearly at 
small qe and reaches a maximum at a certain value of 
qe- This value determines the critical current at which 
the superconducting state become unstable towards the 
normal state. Accordingly, a uniform supercurrent re- 
duces the amplitude of the superconducting order pa- 
rameter, which has been extensively discussed in nodal 
superconductors. 0, H 

One can intuitively understand that Cooper pairs with 
a spontaneous center-of-mass momentum would not oc- 
cur, if the quasiparticle velocity is constant around the 
Fermi surface, for s-wave superconductors due to the 
Galilean invariance. However, once the quasiparticle ve- 
locity strongly depends on the angle around the Fermi 
surface due to Fermi surface nesting, the strength of the 
critical current would vary as a function of the angle, 
(j). This anisotropy of the critical current is an essen- 
tial ingredient for an existence of the FFLO state. In the 
case of anisotropic order parameters, such as an extended 
d-wave order parameter, the anisotropic critical current 
occurs even for a circular Fermi surface. However, we 
found that the region of the FFLO state for this case is 
too narrow to be practically realizable. 

In Fig. 13 we show the direction dependence of the crit- 
ical current for a s-wave superconductor with the tight- 
binding electronic dispersion of the square lattice. The 
strength of the critical current strongly depends on its di- 
rection. It attains its maximum, j^^^, when the current 
is perpendicular to the nesting vector. When the cur- 
rent is applied along a;-direction where the critical cur- 
rent has its minimum, j™™, it is possible to develop a 



spontaneous pair-momentum of q^ to extend the super- 
conducting phase, even above j™'". The development of 
the FFLO state leads to an enhancement of the critical 
current, similar to the enhancement of the critical mag- 
netic field in the previously studied FFLO phase. Below 
we carry out a free energy computation to find a possible 
FFLO state and present temperature and current phase 
diagrams. 

Free energy and Phase diagram: Here we adopt 
the standard BCS formalism to compute the free energy, 
and determine the phase diagram as a function of tem- 
perature and external current. The single particle Green 
function G(iti;„,k) in the presence of finite superflow is 
given byQ 

G"^(iLj„,k) = iujn - Vk •qt/2 + CkP3 + A(k)picri, (2) 

where the single particle dispersion on the square lattice, 
^k — ~t[cos{kx) 4- cos(A;j,)] — /^. The superconducting or- 
der parameter, A(k) ~ Ao/(k), where /(k) represents 
the relative momentum dependence of the order param- 
eter. Pi and Ui are Pauli matrices in particle-hole space 
and spin space, respectively. 

The free energy is then given by 

F-pN ^ -Y, ||L/(^k+) + .f{E^-) - 1] 



^14k'AkAk^, 



kk' 



4 EO^E^^, 



liEy.-)) 



x(l-/(i?k'+)-/(i?k'-)) 

k 

ln(l + /(£;k+))(l + /(i?k-)) 
/3 



(3) 



where E^,± = £;^ ± Vk • qt/2 = VCfc + ^t ± ^k • qt/2. 
The total pair-momentuni, qt is given by the sum of the 
external and spontaneous pair- momenta, qt = qe + q^- 
The pairing interaction 14;k' is assumed to have the form 
of^o/(k)/(k'). 

When an external current is applied along the paral- 
lel axes of the square lattice, qe — q^x (or y), we find 
that the minimum of the free energy as a function of 
the external current shifts from q^ = to a finite value 
of qs above a critical value of g^^. The direction of qs 
is perpendicular to the applied current in order to keep 
constant input and output currents along the direction 
of the applied current. The magnitude of qs depends on 
the magnitude of the external current qe, and it is spon- 
taneously arises so as to maximize the regions of electron 
pairing as indicated by the dark bars in Fig. ^ 

The result of the free energy near zero temperature as 
a function of the external current, qe is shown in Fig. |31 
for several values of qs/qe for the case of s-wave super- 
conductors. Note that above a critical value of g^^, the 
solution with a finite qs has lower energy than that of 
the uniform superconducting state. It is also important 
to note that the superconducting state becomes unsta- 
ble towards the normal state when qe exceeds the critical 
value of q^^ , where the normal state has the lowest en- 
ergy. While our result of Fig. |3| is shown for a s-wave 
superconductor /(k) = 1, the qualitative behavior of the 
phase transition to the FFLO state in the presence of 
current does not depend on the detailed nature of the 
pairing symmetry, /(k). 

The temperature-current phase diagrams for s- and 
d2,2_y2-wave superconductors are shown in Fig. 0] The 
solid line is the second order phase boundary separat- 
ing the normal state and the superconducting state. For 
T/Tc < 0.7 (and 0.6 for a d-wave case), the FFLO state 
which is bounded by the dashed and solid lines is sta- 
bilized. The dotted line indicates the phase boundary 
between the normal and superconducting states in the 
absence of the FFLO state, where the critical current 
at T = is denoted by qco- Note that the FFLO 
state enhances the critical current about 1.7 times for 
n/{2t) = 0.05 at T = for a s-wave superconductor. 
The transition between the FFLO and superconducting 
states is first order for both s- and d-wave superconduc- 
tors at low temperatures. The inset shows the magnitude 
of the spontaneous pair-momentum, qs normalized by the 
external pair-momentum qe as a function of qe/qco- We 
have also analyzed the free energy for a p-wave super- 
conductor with /(k) = sinkx ± i sin ky, and found that 
the qualitative features of temperature-current phase di- 
agram resemble those for the s- and d-wave cases shown 
in Fig. H 

Experimental signatures: Here we assume that the 
superconducting gap and supercurrent will be uniform 
across the system. This condition can be satisfied for 
a sample with the thickness of d <^ £,, where ^ is the 
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FIG. 3: (Color online) The free energy near zero tempera- 
ture vs. qe for various Qs/qe- qe is normalized by qco which 
is the critical current in the absence of the FFLO state. The 
solid, dashed, dotted, and dot-dashed lines are for qs/qe ~ 0, 
0.4, 0.8 and 1.0, respectively. Below qe{= 0.82gco), the uni- 
form superconducting state with Qs = has the lowest energy, 
while above qe{~ 1.6935gco) the normal state has the lowest 
energy. The FFLO states with finite values of Qs are stabilized 
in the window of Qe^ < qe < Qe^ ■ The inset shows the free en- 
ergy vs. qs/qe for various qe near the transition between the 
uniform superconducting and the FFLO states. The solid, 
dashed, dotted lines are for qe/qco = 0.8197, 0.8852, and 
0.9836, respectively. 




FIG. 4: (Color online) A temperature-current phase diagram 
for (a) s-wave and (b) d^2_y2 superconductors. The solid line 
is the second order phase boundary between the normal and 
FFLO states, while the dashed line is the phase boundary 
between the uniform superconducting and FFLO states. The 
dotted line indicates the phase boundary between the normal 
and superconducting states in the absence of the FFLO state, 
where the critical current is denoted by qco a,t T — 0. Tem- 
perature and external current are normalized by the critical 
temperature of Tc and qco , respectively. The inset shows the 
magnitude of the spontaneous pair-momentum qs normalized 
by qe vs. qe/qco 




FIG. 5: (Color online) A schematic diagram for an exper- 
imental set-up. When an external current is applied along 
the axial direction of the tube, where the thickness d <C 5, 
a spontaneous current will be generated along the azimuthal 
direction in the FFLO state. $ denotes the flux induced by 
the spontaneous current. 



coherence length. Below we offer a possible experimental 
set-up to detect the current driven FFLO state. 

Consider a superconductor of cylindrical geometry as 
shown in Fig. [3 If we apply the current along the axial 
direction, the spontaneous current, js will be generated 
along the azimuthal direction in the FFLO state as shown 
in Fig. one of the two- fold degenerate directions of the 
spontaneous current will be chosen. The induced current 
will generate a spontaneous flux, $. We find that a unit 
of flux can be generated when the circumference of the 
tube is approximately 50 nm (ISOxlattice spacing) for 
s-wave superconductors near half-filling with the tight- 
binding electronic dispersion when qs/qe ~ 1- 

The boundary conditions of an order parameter for fi- 
nite size samples deserve some discussion. We suggested 
the cylinder shape of the sample in the above experi- 
ment, so that an order parameter with a finite current, 
A(R) oc e*'i='^ can be stabilized, assuming that we can 
neglect finite size effects on the length of the tube. De- 
pending on boundary, a linear combination of two plane 
waves, such as A(R) ex cos (qs ■ R), can be a possible so- 
lution, which leads to a spatial modulation of the order 
parameter magnitude with a periodicity of 2-k / q^ . 

Recently, the FFLO phase with unequal numbers of 
two hyperfine states of fermionic atoms has been pro- 
posed in the context of cold fermionic atoms. The ex- 
periments have been carried out using Li atoms across a 
Feshbach resonance. The experiment showed that paired 
and unpaired fermions phase separate when the imbal- 
ance is large Isf, but further experiments are needed to 
clarify the nature of the state for a small imbalance of 
densitiesllOj. We suggest that the current driven FFLO 
state can be also realized in optical lattices of cold atoms, 
where tight binding dispersions have been found. 

Summary: We studied the possible existence of the 
FFLO state in the presence of an external current within 



BCS theory. We found that the FFLO state with a spon- 
taneous Cooper pair-momentum can be stabilized in the 
presence of finite currents in superconductors with Fermi 
surface nesting. The anisotropy of the critical current 
due to the nesting plays a crucial role for the existence of 
the FFLO state. A spontaneous pair-momentum is gen- 
erated to extend the superconducting phase and leads 
to an enhancement of the critical current for certain di- 
rections. The FFLO state can be found in anisotropic 
superconductors with a circular Fermi surface, but the 
region of the FFLO state is too narrow to be practi- 
cally realizable. We suggested a possible experimental 
set-up to detect a spontaneous flux induced by a spon- 
taneous current, which will be a direct signature of the 
FFLO state. We also discussed a possible realization of 
the FFLO state induced by a finite current in optical 
lattices of cold fermionic atoms. 
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